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We report on the fluidization dynamics of an attractive gel composed of non-Brownian particles made of 
fused silica colloids. Extensive rheology coupled to ultrasonic velocimetry allows us to characterize the global 
stress response together with the local dynamics of the gel during shear startup experiments. In practice, after 
being rejuvenated by a preshear, the gel is left to age during a time U before being submitted to a constant shear 
rate y. We investigate in detail the effects of both t„, and y on the fluidization dynamics and build a detailed 
state diagram of the gel response to shear startup flows. The gel may either display transient shear banding 
towards complete fluidization, or steady-state shear banding. In the former case, we unravel that the progressive 
fluidization occurs by successive steps that appear as peaks on the global stress relaxation signal. Flow imaging 
reveals that the shear band grows up to complete fluidization of the material by sudden avalanche-like events 
which are distributed heterogeneously along the vorticity direction and correlated to large peaks in the slip 
velocity at the moving wall. These features are robust over a wide range of values of t„ and y, although the 
very details of the fluidization scenario vary with y. Finally, the critical shear rate y* that separates steady- 
state shear-banding from steady-state homogeneous flow depends on the width on the shear cell and exhibits a 
nonlinear dependence with f„,. Our work brings about valuable experimental data on transient flows of attractive 
dispersions, highlighting the subtle interplay between shear, wall slip and aging which modeling constitutes a 
major challenge that has not been met yet. 


I. INTRODUCTION 

Yield stress materials denote assemblies of mesoscopic 
constituents sucb as colloids, droplets or bubbles that display 
macroscopic properties intermediate between those of fluid 
and solid lUEl. At rest, or under low stresses, yield stress 
materials display a solid-like behaviour whereas they flow like 
liquids above a certain threshold referred to as the yield stress 
iSHa. Their solid-like behaviour at rest originates either from 
a densely packed microstructure composed of soft objects m 
such as in dense emulsions 17], jammed microgels ISl fTOl . 
etc., or from the existence of weak attractive interactions that 
bind the constituents together and result in the formation of 
a sample-spanning network HD- In the latter case, solid-like 
properties emerge even at low packing fractions. This defines 
a category referred to as attractive gels lfT2l which encom¬ 
passes various systems such as clay suspensions ifllUBl. car¬ 
bon black gels ESiini and colloid-polymer mixtures IM- All 
sorts of colloids and larger particles of different shapes and at¬ 
traction potential also comply with this definition IIT9ll20l . 

Over the past ten years, the rheological behaviour of attrac¬ 
tive gels has proved to be by far one of the most challenging to 
understand among non-Newtonian fluids. In particular, their 
mechanical properties at rest are strongly time-dependent: at¬ 
tractive gels display a reversible aging dynamics driven by the 
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weak attractive forces between its constituents Ga and that 
can be reversed by shear. As a result, the gel elastic properties 
display a slow logarithmic or weak power-law increase with 
time GijGa, up to complete demixing of the system, which 
is historically referred to as syneresis ll26l . To make matter 
worse, particles are often denser than the surrounding fluid 
which fosters syneresis and may trigger the collapse of the gel 
ETHil]. Such density mismatch is further suspected to in¬ 
fluence the behaviour of these systems under flow, although 
clear experimental evidence is still lacking 122. Nonethe¬ 
less, it is now well established that the behaviour of attrac¬ 
tive gels under external shear involves heterogeneous flows 
that are highly sensitive to preshear history, boundary con¬ 
ditions, and/or finite size effects. For instance, one can em¬ 
phasize the case of Laponite suspensions which steady-state 
flow properties were shown to be influenced by the nature of 
the boundaries: under external shear, smooth walls lead to 
the complete fluidization of the gel and to linear velocity pro¬ 
files, while rough boundary conditions allow the growth of a 
steady shear band ll33l l34l . Moreover, as an illustration of 
the impact of confinement on flows of attractive gels, one can 
mention the spectacular shear-induced structuration observed 
at moderate shear rates and reported in silica suspensions ll35l . 
gels of anisotropic particles ll36l . attractive emulsions ll37ll and 
carbon black gels Il38lj40l . In such cases, the gel rearranges 
to form a striped pattern of logrolling floes aligned along the 
vorticity direction, which origin and formation mechanism are 
still highly debated ll4T1 . 

Beyond the effects of bounding walls and confinement, at¬ 
tractive interactions alone are also responsible for long-lasting 
transients under external shear. On the one hand, experiments 
performed under constant Il42ti46l and oscillatory B3ll47ll48ll 
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stress reveal that the fluidization process of attractive gels ini¬ 
tially at rest, may take up to tens of hours. Experiments cou¬ 
pled to velocimetry have revealed that such a process is mostly 
activated, as evidenced by the Arrhenius-like dependence of 
the fluidization time with the applied shear stress Il43ll45l . and 
that it is strongly heterogeneous in both the vorticity and the 
flow directions ||48l . On the other hand, attractive gels show 
an overshoot in the stress response to shear startup experi¬ 
ments. Such behaviour, which strongly depends on the pres¬ 
hear history, corresponds to the orientation and subsequent 
rupture of the gel microstructure into clusters imiBigiHsi]. 
Beyond the yield point, attractive gels either display homoge¬ 
neous or shear-banded flows depending on the applied shear 
rate and on the boundary conditions Il52ll . However, only a 
handful of studies have investigated the influence of the sam¬ 
ple age, i.e. the duration separating preshear from the start of 
the experiment, on such shear-rate-induced fluidization sce¬ 
nario ll53l . Recently, Martin and Hu have shown on Laponite 
suspensions Il54l that aged-enough samples tend to exhibit 
long-lasting though transient shear-banding, and that shear¬ 
banding may also be trapped by the rapid aging of the non¬ 
flowing band and become permanent. The latter scenario is re¬ 
markable as it strongly differs from the classic shear-banding 
scheme which relies on the mechanical instability of the sam¬ 
ple under scrutiny ll55l 13^ . This study highlights the inter¬ 
play between the sample age and the shear and strongly urges 
to investigate the impact of sample age on the shear-induced 
fluidization scenario in other attractive gels. 

To summarize, attractive gels neither entirely behave as 
solids nor as fluids over a wide range of timescales. In the 
landscape of non-Newtonian fluids, they define a rather singu¬ 
lar category of materials rightfully referred to as the “Twilight 
zone” in the early classification established by L. Bilmes ll57l . 
that was recently adapted to complex fluids by G.H. McKin¬ 
ley 1581. We barely start to understand the (transient) rheology 
of attractive gels and more experiments are obviously needed 
to shed some light on the heterogeneous flows of such highly 
time-dependent materials. 

In the present manuscript we report spatially resolved data 
on the fluidization dynamics of an attractive gel composed 
of weakly attractive non-Brownian particles. Velocimetry 
performed in a concentric cylinder geometry simultaneously 
to shear startup experiments reveals that the steady-state be¬ 
haviour is a subtle function of both the time during which 
the system was left to age before the beginning of the exper¬ 
iment, and of the value j of the applied shear rate. Exten¬ 
sive experiments allow us to build a state diagram of steady- 
state flows in the (j, N) plane. Two distinct regions roughly 
emerge; (/) homogeneous flows for shear rates larger than a 
critical value f* that weakly increases with the aging time, 
and (ii) steady-state shear banding elsewhere. As a key re¬ 
sult of this work, the complete fluidization observed in the 
upper region of the state diagram involves a transient shear 
band that is progressively eroded through a series of dramatic 
fluidization events. These avalanche-like processes show as 
large peaks in the temporal evolution of both the global shear 
stress and the slip velocity measured at the moving wall of 
the shear cell. As further confirmed by two-dimensional ul¬ 


trasonic imaging, this fluidization process is spatially hetero¬ 
geneous and may occur at different locations along the vor¬ 
ticity direction. Einally, for a range of parameters (y, t„) in 
the vicinity of the boundary between the two main regions of 
the state diagram, we observe large fluctuations in the stress 
and slip velocity signals, although the system does not reach 
complete fluidization. Such avalanche-like events are strongly 
coupled to variations in both the width of the shear band and 
the slip velocity. The present work offers to our knowledge 
among the first experimental evidence of local avalanche-like 
fluidization events in a weak attractive gel under shear. It also 
provides an extensive data set to test the relevance of the flow 
stability criteria for shear banding ll59l and stands as a new 
challenge for spatially resolved models ll60l - l^ . 


II. MATERIALS AND METHODS 
A. Experimental setup 

The experimental setup allows us to perform time-resolved 
velocimetry simultaneously to standard rheometry. The rhe¬ 
ological measurements are performed in a polished Taylor- 
Couette (concentric cylinder) cell made of Plexiglas (gap e = 
2 mm) in which the inner Mooney-Couette cylinder of angle 
2.3°, height 60 mm and radius 23 mm is connected to a stress- 
controlled rheometer (ARG2, TA Instruments) and positioned 
at 100 /jm from the bottom of the outer cylinder. A solvent 
trap located at the top of the rotor and a homemade lid are used 
to prevent evaporation efficiently up to about 9 hours. The 
Taylor-Couette cell is embedded in a water tank connected to 
a thermal bath which allows us to keep the sample at a con¬ 
stant temperature T - 23.5 + 0.1°C. 

Local velocity profiles across the gap of the Taylor-Couette 
cell are recorded simultaneously to the global rheology by 
means of two different ultrasonic probes immersed in the wa¬ 
ter tank, which also ensures acoustic coupling to the shear cell. 
The first ultrasonic probe is a single high-frequency focused 
transducer that allows us to record the azimuthal velocity v as 
a function of the radial position r across the gap at the middle 
height of the shear cell i.e. at 30 mm from the bottom. Lull 
technical details on this one-dimensional ultrasonic velocime¬ 
try (ID-USV) technique can be found in a previous publica¬ 
tion ll^ . The second ultrasonic probe consists of a linear ar¬ 
ray of 128 transducers placed vertically at about 15 mm from 
the cell bottom. This transducer array is 32 mm high and gives 
access to images of the azimuthal velocity as a function of the 
radial position r and vertical position z over about 50% of the 
cell height. This two-dimensional ultrasonic velocimetry (2D- 
USV) technique is thoroughly described in ref. ||2l. Both de¬ 
vices can be used simultaneously and roughly face each other 
in the water tank, i.e. they are separated by an azimuthal angle 
of about 180°. While the ID-USV setup has the advantage of 
a better spatial resolution (about 40 //m against 100 //m), only 
the 2D-USV setup allows us to detect and monitor the pres¬ 
ence of flow heterogeneities along the vorticity direction. 

Both velocimetry techniques require that the ultrasonic 
beam crossing the gap of the shear cell is backscattered ei- 
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FIG. 1. SEM images of (a) the stock colloidal suspension (Ludox TM-40) used to prepare the sample, and (b) the non-Brownian aggregates 
composing the sample which form almost instantaneously after mixing the colloidal suspension with a concentrated NaCl solution. The inset 
in (b) shows a single typical aggregate made of fused silica colloids. Images are obtained as follows: both the suspension and the gel have been 
diluted by a factor 1000 (except for the main figure in (b) where the gel has been diluted lOOx), the first one in deionized water and the second 
one in a 6.8% wt. NaCl which corresponds to the final salt concentration in the gel. The samples are then left to dry on a pin stub before being 
imaged, (c)-(e) Bright-field microscopy images representative of the gel microstructure at three different magnifications. The sample has been 
presheared for 2 min at 500 s“' before being imaged about 10 minutes later in a homemade observation cell made of a glass slide topped with 
a cover slip and sealed with UV glue. 


ther by the fluid microstructure itself or by acoustic tracers 
added during sample preparation when the system is acousti¬ 
cally transparent misa. Here, we shall emphasize that the 
microstructure of the sample further detailed below conve¬ 
niently backscatters ultrasound in the single scattering regime, 
which allows us to monitor the fluid velocity in a fully non- 
invasive way. 


B. Sample preparation and rheological properties 

Our Ludox gels are prepared following the same recipe as 
that used by Mpller et al. in ref. Il65]l . As we shall show be¬ 
low, the system is composed of non-Brownian particles made 
of permanently fused silica colloids. These particles are them¬ 
selves reversibly aggregated in brine due to van der Waals 
forces leading to a space-spanning gel with solid-like prop¬ 
erties at rest. 

A stable suspension of silica colloids (Ludox TM-40, 
Sigma-Aldrich, 40% wt. in silica colloids) of typical diam¬ 
eter 20 nm [see Fig. [^a) for a Scanning Electron Microscopy 
(SEM) image of a dilute sample (Supra 55, VP Zeiss)] and 
of pH = 9.0 + 0.5 is first poured without any further purifi¬ 
cation into a 10% wt. deionized aqueous solution of sodium 
chloride (Merck Millipore) up to a mass ratio Ludox:NaCl of 


6:13, corresponding to a final volume fraction of 7% in sil¬ 
ica colloids and a final pH of 7.6 + 0.2. The mixture, which 
instantaneously becomes white and optically opaque, is then 
shaken intensely by hand for 2 min and left at rest for at least 
15 h before being studied. Such a drastic change in the sam¬ 
ple turbidity strongly suggests the rapid formation of aggre¬ 
gates at the micron-scale. Indeed, direct observations using 
different techniques confirm the existence of a much coarser 
microstructure than the initial nanometric silica colloids. On 
the one hand, SEM images of a dried droplet extracted from a 
fresh sample that has been previously diluted in a NaCl solu¬ 
tion unveil the presence of particles which size ranges from a 
few microns up to a hundred microns [Eig.[TJb)]. On the other 
hand, bright-held microscopy images (ECLIPSE Ti, Nikon) of 
the sample neither altered nor diluted further confirm the ex¬ 
istence of these micron-sized particles [Eig. [TJc)-(e)], which 
are stable in time and robust to repeated shear, as confirmed 
by similar observations performed on samples submitted to 
different shear histories (images not shown). 

To account for the formation of such a large scale mi¬ 
crostructure that to our knowledge has not been reported in 
the literature previously l66l . we propose the following 
scenario. Above pH=7, silica colloids are negatively charged 
and bear silanol (Si-OH) and dissociated silanol groups that 
are poorly hydrated ll67l . In most previous works, NaCl is 
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FIG. 2. (a) Elastic modulus G' vs time t after a preshear at jp = 
500 s“' for 2 min. Inset: same data set in semilogarithmic scales. 
Measurements performed under oscillatory shear stress in the linear 
regime with frequency 1 Hz and stress amplitude 0.05 Pa. (b) Flow 
curve (T vs 7 obtained by first decreasing j continuously from 10^ to 
10“^ s“* in 75 logarithmically spaced points of duration dr = 8 s (t) 
and then increasing y over the same range (a). 


added in a relatively small amount (typically 0.05-0.5 M) 
such that electrostatic repulsion is screened leading to the slow 
reversible aggregation of individual colloids up to the forma¬ 
tion of a colloidal gel Il68l470l . Here, we add a much larger 
amount of salt (1.2 M) to the colloidal suspension, which leads 
to an ion exchange where protons are replaced by sodium 
ions Ell Ha. The loss of hydrogen bonding between the col¬ 
loids and the solvent triggers the fast and irreversible aggre¬ 
gation of the silica colloids through the formation of interpar¬ 
ticle siloxane bonds Enca, resulting in the formation of the 
non-Brownian particles described above. Finally, these non- 
Brownian particles also aggregate reversibly due to van der 
Waals forces and form a space-spanning network, i.e. a gel, 
which mechanical behaviour is studied below. Note that such 
microstructure scatters ultrasound efficiently, allowing us to 
use both ID- and 2D-USV without requiring any seeding of 
the sample with tracer particles. 

The rheological features of the gel are displayed in Fig. 

A strong preshear of jp - 500 s ' applied during 2 min fully 
fluidizes the system, which quickly rebuilds once the preshear 
is stopped and subsequently shows pronounced aging. Indeed, 
small amplitude oscillatory shear reveals that the elastic mod¬ 
ulus G' becomes larger than the viscous modulus G" within 


about 20 s [see Fig. 1 in the ESI] and that G' further increases 
logarithmically over more than 2 h [Fig.|^a)]. Such aging dy¬ 
namics are reproducible for a given preshear protocol and lead 
to the formation of a solid-like gel. The latter shows an elastic 
modulus that is nearly frequency independent and a critical 
yield strain of about 7% that weakly depends on the sample 
age [see Figs. 2 and 3 in the ESI]. 

Such a solid-like behaviour also reflects in the nonlinear 
rheology of the gel. Eig. [^b) shows the flow curve cr vs 
7 measured by sweeping down the shear rate from 10 ^ to 
10“^ s * and back up. The flow curve shows an apparent 
yield stress of a few Pascals and displays a complex non¬ 
monotonic shape together with strong hysteresis. Velocity 
profiles recorded simultaneously to the up and down shear 
rate sweeps shown in Eig. [^reveal that wall slip and heteroge¬ 
neous flows are involved over large ranges of shear rates, be¬ 
low 20 s * during the decreasing ramp and up to about 200 s * 
during the increasing ramp [see Eig. 4 in the ESI]. In partic¬ 
ular, the large stress peak observed in Eig. [^b) at 7 ^ 7 s“* 
while ramping up the shear rate is the signature of the partial 
fluidization of the sample which moves as a plug and totally 
slips at both walls for 7 < 7 < 100 s ', before Anally starting 
to flow homogeneously on the reversible branch at high shear 
rates i.e. above about 200 s“*. 


Eurthermore, in Eig. |^b), the shear stress cr shows a min¬ 
imum in between 7 ^ 20 s“* and 40 s“* on the decreas¬ 
ing shear rate sweep, which hints at the existence of a criti¬ 
cal shear-rate jc below which no homogeneous flow is pos¬ 
sible Eaini- The latter result is confirmed by performing 
steady-state measurements at constant applied shear rate start¬ 
ing from the liquid state, i.e. on a fully fluidized sample, in 
order to avoid long-lasting transients that go along with shear 
startup experiments on a sample at rest, and which are at the 
core of section III Here, discrete shear rates of decreasing 
values ranging from 500 s“* down to 0.1 s"' are successively 
applied for at least 300 s each. The flow remains homoge¬ 
neous down to jc ~ 35 s ' below which the sample exhibits 
an arrested band [Eig.j^a)]. Such a value of jc is comparable 
to the one reported in a previous work on very similar Lu- 
dox gels li65l . Moreover, as the shear rate is decreased below 
jc, the size of the fluidized band d decreases roughly linearly 
with 7 [Eig. [^b)] in agreement with the classical “lever rule” 
165117611771. The deviation of d from a straight line could be 
related to the wall slip that is present at the rotor. 


To conclude this subsection, complex cycles of hysteresis 
as the one reported in Eig. |^have been also reported for nu¬ 
merous other attractive gels including carbon black gels na 
and clay suspensions Ii34ll78l . Although significant progress 
has been made in extracting quantitative information from 
hysteresis loops ED, we rather chose to focus our study on 
shear startup experiments in order to fully decouple the fluid 
dynamics from any time-dependent effect related to the ex¬ 
perimental protocol. The present experiments are thus all per¬ 
formed on a sample prepared in a solid-like state, so as to 
investigate the shear-induced fluidization scenario of the non- 
Brownian particle gel. In practice, after preparing the gel in a 
well-defined and reproducible initial state thanks to preshear¬ 
ing, we monitor the stress response of the material to a con- 
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FIG. 3. (a) Steady-state velocity profiles across the gap at different 
shear rates: 7 = 40, 30, 20, 10, 5 and 1 s“' from top to bottom. 
The parameter r denotes the distance to the rotor. These data have 
been recorded starting from large shear rates (500 s“*) and decreasing 
7 by successive steps of long enough duration to achieve a steady 
state at every imposed shear rate, (b) Size of the fluidized band 6 
normalized by the gap width e vs the applied shear rate 7 . The dashed 
line corresponds to the standard “lever rule” with = 35 s“'. 


Slant shear rate over long durations. The results are discussed 
in section Uni 


C. Experimental protocol 

Prior to any shear startup experiment, the sample is pres¬ 
heared at jp - 500 s ' for 2 min in order to erase any previous 
shear history lISOl . We check systematically that the velocity 
profiles during the preshear step become homogeneous across 
the gap within a duration shorter than 2 min. The system is 
then left to rebuild for a duration t„ that ranges between 1 and 
100 min, and during which we apply small amplitude oscil¬ 
latory shear stress (stress amplitude cr = 0.05 Pa, frequency 
/ = 1 Hz) to monitor the evolution of the gel linear viscoelas¬ 
tic properties and make sure that for a given value of the 
gel reaches an initial state that is reproducible from one exper¬ 
iment to another. Finally, a constant shear rate 7 is applied to 
the material over a long duration (> 1000 s) while we monitor 
the stress response together with local velocity profiles. 


III. RESULTS 

A. Global rheology and one-dimensional velocity profiles 

We first fix the duration between the preshear and the start 
of the experiment to t„ = 60 min and discuss the material 
response for shear startup experiments performed at different 
shear rates. For shear rates lower than about 150 s ', the stress 
relaxes smoothly towards a constant value as illustrated in fig¬ 
ure a) for an experiment performed at 7 = 100 s '. Veloc¬ 
ity profiles acquired simultaneously reveal that the material 
remains at rest near the stator, while it flows in the vicinity 
of the rotor where it shows strong wall slip [Fig. [^b)]: the 
sample displays steady-state shear banding. For larger shear 
rates, one observes a completely different scenario, both from 
global and local measurements. The stress displays a series 
of spike-like relaxation events, during which stress increases 
of small amplitude are followed by large drops [Fig. Qa)]. 
Local measurements show that the gel flows heterogeneously. 
However, the shear band is only transient and the system al¬ 
ways ends up being homogeneously sheared at steady state 
[Fig.gc) -(e)]. Qualitatively, the width of the shear band in¬ 
creases while strong wall slip is observed at the rotor, until 
the whole gap is fluidized and wall slip becomes negligible. 
This scenario is robustly observed at different shear rates and 
full fluidization occurs sooner for larger shear rates [Fig. Qa)] . 
A quantitative analysis of the velocity profiles is further pro¬ 
posed in section [III C| 

To evidence the impact of the aging time t„ on the material 
response. Figs. Qf)-(j) show similar experiments performed 
for different t„ at the same shear rate j - 100 s"'. For long 
waiting times between the preshear and the start of the experi¬ 
ments, e.g. t„ - 100 min, one observes a smooth stress relax¬ 
ation here again associated with steady shear-banded velocity 
profiles [Fig. j^f) and (g)]. Decreasing t„ to 30 or 5 minutes 
leads to a gel of lower elastic modulus. Local measurements 
further reveal that these weaker gels go through a transient 
shear-banding regime and that in both cases the steady state 
is homogeneous [Fig. Qh) and (i)]. Here, unlike the case of 
transient shear banding reported in Fig. j^a) for large shear 
rates in samples left to age for t„ - 60 min where fluidiza¬ 
tion corresponds to a long series of successive stress relax¬ 
ations, fluidization proceeds in a single stress drop together 
with small noisy fluctuations [Fig.j^f)]. Finally, one observes 
that very young gels {t„ - 1 min) barely show any heteroge¬ 
neous velocity profile during startup flow and reach a homoge¬ 
neous steady state within a few tens of seconds [Fig.|^)]. In 
summary, the longer the sample ages after preshear, the more 
likely it is to exhibit a long-lasting heterogeneous fluidization 
scenario or to display steady shear banding. 


B. Flow state diagram 

To go one step further and get an overall picture of the ma¬ 
terial response, we have performed extensive shear startup ex¬ 
periments by varying systematically both the aging time (1 < 
t„ < 100 min) and the applied shear rate (10 < 7 < 250 s“'). 
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FIG. 4. (a) Shear stress response cr(f) for different applied shear rates [color, y (s“')]: [— , 100]; [— , 150]; [— , 200]; [—, 250]. (h)-(e) Velocity 
profile v(r) across the gap, where r is the distance to the rotor at different times [symbol, time (s)]: [O. 200]; [□, 400]; [A, 800]; [V, 1400]. 
Each color corresponds to each applied shear rate in (a). The rotor velocity corresponds to the upper bound of the vertical axis. The sample is 
left to age during t„, = 60 min before each experiment, (f) Shear stress response cr(f) for different waiting times t„ [color, t„ (min)]: [— , 100]; 
[— , 30]; [— , 5]; [—, 1]. Experiments performed at y = 100 s“*. (g)-(j) Velocity profile at different times [symbol, time (s)]: [O, 30]; [□, 400]; 
[A, 800]; [V, 1100]. Each color corresponds to each value of t„ in (f). The rotor velocity corresponds to the upper bound of the vertical axis. 


The entire data set is summarized in the flow state diagram 
pictured in Fig. where the different symbols code for the 
three different scenarios that we may distinguish at the end of 
each startup experiment lasting typically 1500 to 3600 s. In 
steady state, the attractive gel may either fully fluidize after 
a transient phase possibly involving strong fluctuations (♦ in 
Fig.ig, as observed in the upper part of the state diagram for 
large shear rates quite independently of the sample age, or ex¬ 
hibit shear banding. The boundary between these two regimes 
defines a critical shear rate f* that depends in a nonlinear fash¬ 
ion on the sample age t„. Furthermore, we can discriminate 
between two types of behaviour when shear banding is ob¬ 
served. 

First, steady shear banding (• in Fig.|^, for which both the 
width 6 of the shear band and the slip velocity Vj at the ro¬ 
tor display negligible fluctuations in steady state, is observed 
in particular at low aging times t„ < 20 min and under low 
enough shear rates, typically below 100 s '. In that case, the 
width of the arrested band may decrease (f„ = 1 s) or remain 


constant {t„ = 10 s and 100 s) for increasing shear rates. 

Second, we have also observed unsteady shear banding 
(♦ in Fig. 1^. In this case both global and local measure¬ 
ments display significant fluctuations in steady state. These 
fluctuations are strikingly similar to those observed during the 
transients leading to complete fluidization in the upper part of 
the diagram. However here the material never gets entirely 
fluidized and the shear band width 6 does not show a system¬ 
atic evolution towards 6 - e so that an unsteady shear band 
persists in steady state at least within the finite duration of the 
experiments. These fluctuations are investigated in more de¬ 
tails in the next section. 


C. Linking global rheology to the local dynamics 

In this section, we focus on the strong fluctuations that are 
observed (i) during transient regimes leading to complete flu¬ 
idization and (ii) during unsteady shear-banded flows at steady 
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FIG. 5. Flow state diagram in the (sample age f„,, applied shear rate 
y) plane. In steady state, the gel may either be fully fluidized (♦, 
blue region) or display shear banding. Steady shear banding is rep¬ 
resented in (•) symbol which size encodes the portion of the gap that 
is being sheared. Unsteady banding, which denotes flows where the 
band width and the slip velocity at the rotor display significant fluc¬ 
tuations, is represented by (♦). 


state. Global rheological data together with the detailed anal¬ 
ysis of the corresponding local measurements are reported in 
Figs. HI and |7] respectively (see also corresponding Movies 1 
and 2 in the ESI). Local data are analyzed as follows: lin¬ 
ear fits of the velocity profiles in the fluidized part of the gap 
are used to estimate the local shear rate and to extrapolate the 
sample velocity v(0) at the rotor. The width 6 of the shear band 
is obtained as the abscissa of the intersection between the fit 
and the zero velocity axis, while the slip velocity is given by 
Vi = vo-v(O), where vq is the rotor velocity. Finally, the global 
shear rate is defined as v(0)/e. 

We first discuss Fig. which shows a case of full fluidiza¬ 
tion for j = 150 s“' on a gel left to age for t„ = 30 min. 
The stress relaxes by successive jumps up to full fluidization 
which occurs at f = t/ ^ 507 s and after which cr remains 
roughly constant [Fig. |^a)]. The temporal evolutions of d 
[Fig. I^b)] and of the local shear rate [Fig. |^c)] show that 
fluidization occurs by successive spatial “avalanches” that are 
directly correlated to the stress drops. During the whole flu¬ 
idization process, the slip velocity at the rotor V j decreases to¬ 
wards a negligible value that is reached at f = t/, diminishing 
by jumps that are in phase with the stress evolution. 

Figure |7] focuses on spatiotemporal fluctuations observed 
during steady-state shear-banded flows for the same aging 
time (t„ = 30 min) as in Fig. but sheared at a lower shear 
rate (j = 50 s“'). The stress exhibits periods of slow increase 
followed by rapid drops [Fig.j^a)]. Within an hour, about half 
the gap gets fluidized and the size of the fluidized band shows 
pronounced increase during short period of times, that are syn- 



FIG. 6. Full fluidization observed for f = 150 s“' and 

G = 30 min. (a) Shear stress cr (—) and slip velocity 

Vj (— ) vs time. (b) Width S of the fluidized shear band 
normalized by the gap width e vs time. (c) Local shear 
rate within the shear band (—) and global shear rate (— ) 
vs time. The horizontal dotted line indicates the shear rate applied 
by the rheometer. In (a)-(c), the vertical dashed line indicates the 
fluidization time t/. 


chronized with the stress drops [Fig. |^b)]. The dynamics of 
the fluid at the rotor is strongly correlated to the global fluctu¬ 
ations, as evidenced by the sudden peaks of the slip velocity 
when the stress drops. Note that, although the slip velocity 
decreases in average over the whole duration of the experi¬ 
ment, it remains at a high level of about 20%, in noticeable 
contrast with the fully fluidized scenario described in the pre¬ 
vious paragraph. Such oscillatory dynamics in the vicinity of 
the rotor are reminiscent of stick-slip. Indeed, the fluidized 
band shows time interval during which the local shear rate re¬ 
mains constant [Fig. |^c)], while the band width slowly de¬ 
creases. As a result, the stress slowly builds up, until the 
fluidized band experiences a large slip event at the moving 
wall and gets rejuvenated. These dynamics look similar to the 
stick-slip motion reported in clays in the pioneering work of 
Pignon et al. IISTII . However, in the latter case, stick-slip oc¬ 
curs along fracture planes located in the bulk sample, while 
here slip at the wall appears to play a key role. 

Finally, we shall emphasize that in both types of scenar¬ 
ios, the peaks in both the stress and the slip velocity are only 
seen in the presence of shear banding. Although it is not clear 
which of the shear banding or of the slip at the rotor is the 
cause of the oscillations, these peculiar dynamics hint toward 
a subtle flow-microstructure coupling that certainly deserves 
more investigation. 
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FIG. 7. Unsteady shear banding observed for y = 50 s“' and 
tn- = 30 min. (a) Shear stress cr (—) and slip velocity Vj (— ) vs time, 
(b) Width S of the fluidized shear band normalized by the gap width 
e vs time, (c) Local shear rate within the shear band (—) and global 
shear rate (— ) vs time. The horizontal dotted line indicates the shear 
rate applied by the rheometer. 



FIG. 8. (a) Stress cr vs time during a single stress drop event extracted 
from shear startup experiment during which the material is fully flu¬ 
idized (y = 150 s“*, t„ = 10 min). Inset: stress vs time for the whole 
experiment. The signal pictured in the main graph appears in red. 
(b) Spatiotemporal diagram of the velocity data v(r, t) as a function 
of position r and time t. The radial position r is measured from the 
rotating inner wall. Data obtained with ID-USV. (c) Spatiotemporal 
diagram of the velocity data v(ro, z, t) as a function of the vertical po¬ 
sition z and time f at tq = 0.2 mm. Data obtained with 2D-USV. The 
vertical position z is measured from the top of the transducer array. 
In both (b) and (c), the fluid velocity is color coded in mm.s^'. 


D. Local scenario within an avalanche 

To get better insight on the avalanche-like fluidization sce¬ 
nario, we now focus on the single stress drop shown in 
Fig. I^a). This event is extracted from a shear startup ex¬ 
periments performed at y = 150 s ' for t^, - 10 min and 
the steady-state flow corresponds to a fully fluidized sample 
[inset of Fig. a)]. The velocity profiles recorded simultane¬ 
ously with ID-USV are plotted as a spatiotemporal diagram 
in Fig.[^b). Quite surprisingly the one-dimensional flow pro¬ 
files before and after the stress drop are very similar as they 
all show shear localization over about half the gap. By the 
time the stress reaches its maximum value, about 80% of the 
gap is sheared. Then shear abruptly localizes again over about 
1 mm close to the rotor at the beginning of the stress relax¬ 
ation. Therefore, although the strong fluctuations of v(r, t) ob¬ 
served before the stress peak appear to be correlated to those 
of cr(f), the drop of about 15% in the stress value cannot be 
explained by these one-dimensional data in a straightforward 
manner. We emphasize that the ID-USV measurements are 
performed at a given height of the Taylor-Couette cell so that 
they may reflect the evolution of global rheology only if the 
flow is homogeneous along the vorticity direction. As a matter 
of fact, velocity profiles v(r, z, t) recorded simultaneously over 
the whole height of the Taylor-Couette cell through 2D-USV 
demonstrate that the flow is strongly heterogeneous in the ver¬ 
tical direction z. Figure [^c) shows a spatiotemporal plot of 
the fluid velocity along the vertical axis z at a fixed radial po¬ 


sition ro - 0.2 mm close to the rotor (see also Movie 3 in the 
ESI). Depending on the position along the z-axis, the material 
can simultaneously be either solid-like as evidenced by areas 
of very low velocities in Fig. [^c) (see for instance at the top 
of the region of interest for z < 5 mm) or fluid-like and flow at 
high velocities as observed for z - 8 and z - 30 mm. Note that 
the ID-USV measurements are performed at zo - 15 mm and 
are fully consistent with the 2D-USV data. Our measurements 
further reveal that the particular avalanche-like event analyzed 
here only corresponds to partial fluidization: at the end of the 
stress drop, the flow is still heterogeneous, even though ar¬ 
rested regions for t < 870 s has given way to flowing regions 
arranged in a vertically banded structure for f > 890 s. In fact, 
the inset of Fig. a) shows that several subsequent stress drop 
events still have to take place before complete fluidization of 
the sample. 

To test the existence of a local fluidization scenario that 
would be generic to all avalanche-like events, we analyze ad¬ 
ditional experiments in the full fluidization region of the state 
diagram at y = 180 and 200 s *. The temporal evolution of 
the shear stress is reported respectively in Figs.|^a) andfTOta) 
and the corresponding 2D-USV measurements v(ro,z, f) are 
displayed as spatiotemporal diagrams in Figs.j^b) and [T0|b) 
for a radial position ro = 0.7 mm (see also Movies 4 and 5 
in the ESI). Here again, the temporal evolution of the veloc¬ 
ity field is strongly heterogeneous along the vertical direction 
in both experiments. Eurthermore, despite these two experi- 


8 









































D Local scenario within an avalanche 
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FIG. 9. (a) Stress cr vs time during a single stress relaxation episode 
extracted from shear startup experiment at the end of which the ma¬ 
terial is fully fluidized (7 = 180 s“', t„ = 1 min), (b) Spatiotemporal 
diagram of the velocity data v(ro, z, t) as a function of the vertical po¬ 
sition z and time f at ro = 0.7 mm. Data obtained with 2D-USV. The 
vertical position z is measured from the top of the transducer array. 
The fluid velocity is color coded in mrn.s^*. 


ments are performed at comparable shear rates, the local flu¬ 
idization scenario is strikingly different. In the shear startup 
experiment reported in Fig. the lower part of the region of 
interest is fully fluidized after the first stress relaxation while 
the upper part of the sample necessitates two supplemental 
avalanches to in turn fully fluidize. On the other hand, the flu¬ 
idization process shown in Fig. 10 starts from the upper part 
of the region of interest before extending to the lower part 
of the Taylor-Couette cell. Here, each avalanche-like event 
involves large pieces of the sample with a typical vertical ex¬ 
tension of about 5 mm [see events marked by white dotted 
lines in Fig.[T0|a)]. Some of these events only show up in the 
stress response as very small peaks while their local signature 
is much more impressive (see, e.g., the event at f ^ 600 s). 
For t ^ 400 to 500 s, the sample even appears to fluidize, or 
at least to be set into motion, through regular steps occuring 
from top to bottom with a characteristic time of 10-20 s. The 
experiments shown in Figs.[^and[T^allow us to conclude that 
the peaks in the stress signal may encompass very different 
local scenarios. Since both experiments were performed with 
the same aging time t„ - \ min, these results illustrate the 
high sensitivity to “initial conditions,” i.e. to the possibly dif¬ 
ferent arrangement of the heterogeneous microstructure after 
the aging process, and the subtle interplay between aging and 
shear banding with no systematic failure scenario along the 
cell height. 



t (s) 


FIG. 10. (a) Stress cr vs time during a single stress relaxation episode 
extracted from shear startup experiment at the end of which the ma¬ 
terial is fully fluidized (7 = 200 s“', t„ = 1 min), (b) Spatiotemporal 
diagram of the velocity data v(ro, z, t) as a function of the vertical po¬ 
sition z and time t, at ro = 0.7 mm. Data obtained with 2D-USV. The 
vertical position z is measured from the top of the transducer array. 
The fluid velocity is color coded in mm.s^'. 


IV. DISCUSSION 

On the one hand, we conclude from section |IIB| that the 
gel under scrutiny displays a discontinuous yielding transi¬ 
tion that is very similar to the one reported in other attractive 
gels 165] im [821. Indeed, one observes that, when the ini¬ 
tial condition is a completely fluidized state, applying a se¬ 
ries of shear rates of decreasing values leads to the growth 
of a steady-state arrested band below a critical shear rate jc- 
The value of 7 ^ is inherent to the sample and results from 
the flow instability at low shear rates. This is fully consis¬ 
tent with previous results on Ludox gels, li^ On the other 
hand, the shear startup experiments reported in section jUIj and 
performed on the sample prepared in the solid state allow us 
to identify a second critical shear rate 7 * > jc- For applied 
shear rates such as 7 ^ < 7 < 7 *, the sample gets only partially 
fluidized and exhibits steady-state shear banding, while long- 
lasting yet transient shear banding is observed for 7 > 7 *. 
Such a behaviour can be interpreted in the framework recently 
described by Martin and Hu on Laponite samples ll54ll . Un¬ 
der an external shear, the sample in the vicinity of the rotor 
that is in initially at rest gets fluidized, while the sample close 
to the stator remains at rest and keeps aging. Depending on 
the intensity of the shear rate, such heterogeneous velocity 
profile may either be trapped by the sample aging and become 
permanent as observed below 7 *, or become homogeneous af¬ 
ter a transient shear-banding regime which duration decreases 
for increasing values of 7 . Such an argument allows us to un¬ 
derstand that 7 * will be sensitive to the gap size and to the 
cell geometry, which we have verified by performing identi- 
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cal shear startup experiments with different gap widths. These 
experiments show that f* decreases as the gap size is reduced 
[see Fig. 6 in the ESI]. Therefore, steady-state shear banding 
trapped by aging should be distinguished from shear banding 
resulting from intrinsic flow instability which should not de¬ 
pend on the geometry. Finally, we have shown in section [IIlB| 
that f* is a nonlinear function of the sample age t„, again in 
contrast with %. 

Regarding the transient fluidization dynamics, our work un¬ 
ravels the existence of very characteristic peaks in the global 
rheological data. These peaks correspond to local avalanches 
associated with the abrupt fluidization of shear-banded ve¬ 
locity profiles. An avalanche proceeds in two steps. First, 
the sample ages as evidenced by the slow increase of the 
stress indicative of progressive consolidation, while the shear 
band remains roughly fixed. Second, the sample suddenly 
fluidizes before localizing again, at least partially, while the 
stress drops and the shear rate increases. This scenario strik¬ 
ingly recalls the transient shear banding reported in Laponite 
clay suspensions 1541 and fits well with the stability criterion 
recently proposed by Fielding et al. Il59l |62l . In practice 
though, the physical origin of the avalanche-like and succes¬ 
sive stress relaxation events remains to be determined. Aging 
is more pronounced in our system than in the attractive gels 
that have been studied previously ll22] (831 (841 as evidenced 
from the large values of dG'/dlogf after preshear [see inset 
of Fig. I^a)]. Therefore, syneresis driven by the aggregation 
and/or sedimentation of the colloidal floes due to their density 
mismatch with the surrounding fluid could also play a role and 
interfere with the traditional picture of a competition between 
aging and shear rejuvenation ll85l . As syneresis is negligible 
in most of the colloidal gels that have been the topic of rheo- 
physical studies so far, it could also explain why such stress 
oscillations have, to our knowledge, never been reported in 
the literature. 

An alternative explanation for the stress oscillations could 
be related to confinement. Indeed, as the size distribution of 
the fused silica aggregates is wide and extends up to 100 yum, 
the sample mesostructure most likely involves aggregates 
which size becomes comparable to that of the gap, at least 
for long enough rest durations In this framework, stress 
oscillations would result from a competition between shear- 
induced structuration as described in the introduction and the 
strong aging of the sample. Such an interpretation would also 
account for the stick-slip like motion of the fluidized band at 
the moving wall. Nonetheless, despite systematic monitoring 
of the sample during shear startup experiments, no structura¬ 
tion or spatial pattern could be observed. Moreover, supple¬ 
mental shear startup experiments in narrower gaps show that 
for a given shear rate, stronger confinement leads to the dis¬ 
appearance of the stress oscillations and to homogeneous ve¬ 
locity profiles [see Fig. 6 in the ESI]. This last result strongly 
suggests that confinement alone cannot account for stress os¬ 
cillations. 

Finally, the present study has focused on experiments per¬ 
formed under smooth boundary conditions, revealing the pres¬ 
ence of strong wall slip associated with heterogeneous, shear- 
banded flows while fully fluidized states show negligible wall 


slip. Yet the roughness and/or chemical nature of the walls 
are known to have a crucial influence not only on rheological 
measurements but also on local flow both close to the walls 
and in the bulk Il3^ 1^ l45l 1^ 1^. Therefore the influence 
of boundary conditions on the complex fluidization scenario 
reported here appears as the next key question to address in 
future work. As a first step, we report preliminary tests on the 
role of boundary conditions on the above results in the ESI. 
The flow curve measured under “rough” boundary conditions 
in a sand-blasted Plexiglas Taylor-Couette cell with typical 
roughness of 1 yum (to be compared to a few tens of nanome¬ 
ters for the “smooth” cell used so far) shows a smaller, yet 
significant hysteresis [see Fig. 1(b) in the ESI] as well as wall 
slip at low shear rates (see Fig. 5 in the ESI). Furthermore, 
one can see a strong difference between the velocity profiles 
recorded simultaneously to the flow curve respectively with 
the rough and smooth boundary conditions (compare Fig. 4 
and 5 in the ESI) although the surface roughness of the rough 
boundary does not match exactly the size of the microstruc¬ 
ture. In particular the transient fluidization episode reported 
around 7 s“^ with smooth boundary conditions [Fig.j^b) and 
Fig. 4 in the ESI] is absent with rough walls (see Fig. 5 in the 
ESI). These preliminary results illustrate the strong impact of 
the boundary conditions and urge for systematic experiments 
so as to quantify the impact of the boundary conditions espe¬ 
cially on the state diagram reported in Fig.[^ 

V. CONCLUSION 

We have investigated the local scenario associated with the 
shear-induced fluidization of an attractive gel made of non- 
Brownian particles. We have identified a critical shear rate y* 
that separates steady shear-banded flows from full fluidization 
and that exhibits a nonlinear dependence with sample age t„. 
This critical shear rate y* is much larger that the one signaling 
flow instability in experiments starting from fluidized states 
and depends on the shear cell geometry as well as possibly on 
the preshear intensity imposed prior to the shear startup ex¬ 
periment. As a key result, we have shown that for shear rates 
larger than y*, the fluidization of the sample involves succes¬ 
sive local avalanche-like events that are heterogeneously dis¬ 
tributed along the cell height and which dynamics is strongly 
coupled to both the slip behaviour at the wall and the width 
of the shear band. Such avalanches appear in the stress sig¬ 
nal as peaks, which individual properties are strongly remi¬ 
niscent of stick-slip phenomena. Future work will focus on 
the early stage of shear startup experiments, and in particular 
on the stress overshoot that occurs before the stress relaxation, 
as well as on the influence of the boundaries and confinement 
on the fluidization scenario. 
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SUPPLEMENTAL MATERIAL 


Appendix A: Supplemental movies 

Five movies are provided as supplemental materials to 
Figs. 5-9 of the main text. All movies illustrate the tempo¬ 
ral evolution of the global rheology and of the ID and/or 2D 
velocity profiles recorded simultaneously to the rheology dur¬ 
ing successive avalanche-like events. 


Appendix B: Supplemental figures 

Supplemental Fig. [TT]^a) displays the evolution of the elas¬ 
tic and viscous moduli of the sample after preshear. The gel 
rebuilds quickly as evidenced by the fact that the elastic modu¬ 
lus is larger than the viscous modulus about 20 s after preshear 
has been stopped. 

Supplemental Fig. [TTJb) shows flow curves cr vs f mea¬ 
sured by rapidly sweeping down the shear rate from 10^ to 
10“^ s ' and back up. The flow curves have complex non¬ 
monotonic shapes together with strong hysteresis under both 
smooth and rough boundary conditions. 

Supplemental Fig. [T^ illustrates the result of a frequency 
sweep performed at a constant strain (y - 0.5 %) from 10 



A 
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10 10 10 10 10 
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FIG. 11. (a) Elastic (o) and viscous moduli (A) vs time after a pres¬ 
hear at fp = 500 s“' for 2 min (f - 1 Flz, cr = 0.05 Pa). Same data 
as in Fig. 2(a) in the main text, (b) Flow curve, shear stress cr vs. 
applied shear rate y, obtained by decreasing (V) and increasing (A) 
shear rate in order from 10“^ to 10^ s“' with a waiting time of 8 s per 
point. Black (resp. red) symbols correspond to results in a smooth 
(resp. rough) geometry. 


to 0.02 Hz, and for different aging times = 30, 60 and 
100 min after stopping the preshear. The elastic modulus G' 
is independent of the frequency and increases with the aging 
time tu, in agreement with Supplemental Fig. 11 The viscous 
modulus shows a broad minimum with no sign of increasing 
part at low frequencies, a feature that is typical of soft glassy 
systems m. 

Supplemental Fig. 13 shows the results of a strain sweep 
performed on the gel at a constant frequency / = 1 Hz for 
different aging times after preshear. The gel yields at a 
constant strain of about 7 % that is roughly independent of the 
aging time t„. 

Supplemental Figs. 14 and 15 show the rheological and 
2D-USV data obtained in smooth and rough boundary con¬ 
ditions respectively, by first decreasing the shear rate and then 
increasing it back again. The corresponding flow curves are 
shown in Fig. 1(b) of the present ESI. Here we also show spa- 
tiotemporal plots of the velocity profiles as a function of the 
radial position r at a given height (zo - 15 mm) in the Taylor- 
Couette cell [Supplemental Figs. [I4|^b) and[T5|[b)], and as a 
function of the vertical position z, at a given radial position 
inside the gap (ro = 0.5 mm) [Supplemental Figs. [T4|c) and 
[lie)]. Velocity data are shown only for y <460 and 340 s * 
in smooth and rough boundary conditions respectively, due to 
the limitation in the ultrasonic pulse repetition frequency to 
20 kHz, which sets an upper bound on the measurable veloci¬ 
ties JS). 

For smooth boundary conditions, the flow curve exhibits a 
large hysteresis, mainly localized in a narrow range of shear 
rates [Supplemental Fig. [TTJb)]. During the decreasing ramp 
of y, the sample which is first sheared homogeneously at large 
shear rates, enters a total wall slip regime at f ^ 200 s, i.e., for 
■y ~ 20 s“* [Supplemental Fig. 14 b)]. This regime persists 
until y is increased back up again. The subsequent fluidization 
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FIG. 12. Elastic (A) and viscous (V) moduli vs frequency /. The 
frequency sweep is performed at a constant strain y = 0.5 % and 
at different aging times t„ after preshear, [color, t„ (min)]: (-, 30); 
(-, 60); (-, 100). 
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FIG. 13. Elastic (o) and viscous (A) moduli vs the strain amplitude. 
The strain sweep is performed at a fixed frequency / = 1 Hz with 
a waiting time of 8 s per point, and at different times G after the 
preshear, [color, f„, (min)]: (-, 30); (-, 60); ( , 100). The vertical 
dashed line at y ~ 7 % emphasizes the yield point defined as the 
point at which G' = G". 



and that the local behavior of the sample is roughly homoge- 
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FIG. 15. Same as Supplemental Fig. [T^ for a Taylor-Couette cell 
with rough boundary conditions. 



FIG. 14. Taylor-Couette cell with smooth boundary conditions, (a) 
Shear rate y and shear stress <t vs time t obtained by first decreas¬ 
ing continuously y from 10^ to 10“^ s“* in 75 logarithmically spaced 
points of duration dt = 8 s each, and then increasing y over the same 
range, (b) Spatiotemporal diagram of the velocity data v(r,zo,l) as 
a function of the distance r to the rotor and time ?, at zo - 15 mm. 
(c) Spatiotemporal diagram of the velocity data v(ro,z,t) as a func¬ 
tion of the vertical position z and time f, at ro = 0.5 mm. The verti¬ 
cal position z is measured from the upper boundary of the 2D-USV 
probe. The fluid velocity is color-coded in mm.s^'. 


of the material is abrupt and occurs in a narrow range of shear 
rates, 1 <y< 10 s \ which coincides with the bratal increase 
of the global stress. Supplemental Fig.[T4|b,c) reveal that the 
subsequent drop in the stress corresponds to flow arrest, i.e. 
total wall slip, until the sample fluidizes again for 7 > 100 s“^ 


neous along the whole height of the cell during both ramps. 


Supplemental Fig. 15 illustrates the exact same protocol 


as that reported in Supplemental Fig. but performed in 
a rough Taylor-Couette cell. As also seen in Supplemental 
Fig.[TT]^b), the rheological hysteresis is far less pronounced in 
comparison to the data obtained with smooth boundary con¬ 
ditions. Here the flow remains homogeneous along the ver¬ 
tical axis during both ramps [Supplemental Fig. [T5|c)]. The 
velocity profiles at a given height in the Taylor-Couette cell 
also display very different behavior from the smooth case. 
In the decreasing ramp of y, the stress plateaus in between 
f ~ 100 s ' and y ~ 1 s ', while shear localizes at the rotor. 
Below y ~ 1 s“*, the stress displays a kink towards a sec¬ 
ond plateau, while the sample exhibits a plug flow. Increasing 
y back up, the plug flow persists up to shear rates of a few 
10 s *. Above the latter value, the fluidization of the sample 
proceeds from the rotor and involves transient banding. Full 
fluidization is not available in Supplemental Fig.[T5][b) for the 
technical reason mentioned above. 


Supplemental Fig. 16 shows the influence of the gap size on 
the fluid response to a shear startup experiment performed at 
y = 50 s * on a sample left at rest for = 10 min. In a gap of 
size e = 2 mm, the same value as in the main text, the stress 
response displays a large number a peaks together with un¬ 
steady shear banding over 2400 s [Supplemental Fig.[T 6 tb)]. 
Decreasing the size of the gap changes radically the sample 
behavior. In a gap twice as small (e - 1 mm), the stress re¬ 
sponse exhibits less peaks and the sample fluidizes entirely in 
about t - 200 s [Supplemental Fig. [Thtc)]. Decreasing the 
gap by a factor of 4 (e = 0.5 mm) leads to a smooth stress 
response and to homogeneous velocity profiles from the be¬ 
ginning of the experiments, with significant slip at the moving 
boundary [Supplemental Fig.fThtd)]. 


14 





























B SUPPLEMENTAL FIGURES 




FIG. 16. (a) Stress response cr(f) to a shear startup experiment performed at y = 50 s“' for different gap sizes: [color, gap (mm)] = [— , 2], 
[— , 1], [—, 0.5]. (b)-(d) Velocity profile, where r is the distance to the rotor at different times in (a) [symbol, time (s)]: [O, 100]; [□, 800]; [A, 
1200]; [V, 2000]. The rotor velocity corresponds to the upper bound of the vertical axis. The sample is aged during = 10 min before each 
experiment. 
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